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Abstract 


Fragility curves are derived from fragility function that indicates the probability of damage of structure 
due to earthquake as a function of ground motion parameter. It helps to predict the level of structural 
damage and consequently reduce the seismic risk in specific ground motion. In this scenario, this study 
is focused on the construction of fragility curve of institutional reinforced concrete (RC) building of 
Pokhara University. For this, the building of School of Health and Allied Science (SHAS) is 
considered as a guiding case study. For the numerical analysis, the study building blocks are modelled 
in finite element-based software. The non-linear static and linear dynamic analyses are employed for 
numerical analysis. In dynamic analysis, building models are subjected to the synthetic accelerograms 
of the 2015 Gorkha earthquake. Based on the analyses, the analytical fragility curves are plotted in 
terms of probability of failure at every 0.1 g interval of peak ground acceleration (PGA) with log 
normal distribution. Finally, the results are highlighted for different seismic performance level in 
buildings: slight damage, moderate damage, extensive damage and complete damage for the 
earthquake of 475 years return period. 


Keywords: Non-linear analysis, Fragility function, Seismic performance, Reinforced concrete 
building 


1. Introduction: 


Nepal is located in the Himalayan region which 
has a long history of moderate to strong seismic 
events. The main cause of earthquake in the 
Himalayan belt in due to the collision between 
Indian and Eurasian plate [1, 2]. Pokhara valley is 
located near the main central thrust, that one of the 
most active fault which caused the repeated 
powerful earthquakes on past. Geological 
variation, terrain characteristics, urbanization of 
city area, low construction materials and technics 
as well as age of buildings in Pokhara valley are 
prime causes of structural vulnerability, lead to 
dead of people and loss of economy [3]. 


Initially, the institutional buildings in Nepal were 
constructed with adobe, brick masonry, stone 
masonry, composite masonry, and timber frame 
structures. After the increasing maturity of 
education system as well as construction 
industries in Nepal, RC building construction has 
started from late 1970s [4]. Currently, institutional 
buildings like adobe, brick masonry, stone 
masonry and timber are replaced by reinforced 
concrete construction. In the past, most of the 
buildings face several major earthquakes in the 
year of 1988, 2011 and 2015. After these major 
earthquakes, many institutional buildings were 
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functional without assessing any seismic 
evaluation. This might be a prime concern for 
future catastrophic. 


Researchers highlighted the damages of 
institutional buildings which leads large number 
of casualties and economic loss in their study [5- 
7]. Miyamoto et al., [5] highlighted the damage 
and collapse of school building in Indonesia due 
to Sumatra earthquake. Augenti et al., [6] 
discussed damage and collapse of school building 
in Italy. In china, Sichuan earthquake damaged 
the thousands of buildings and injured millions of 
students [5]. More than 6,000 school building 
were damaged that lead 17,000 students and 900 
teachers died by 2005 and 2008 earthquake in 
Pakistan [7]. The 2015 Gorkha earthquake in 
Nepal damaged more than 8200 institutional 
buildings and an estimated loss of nearly 285 
million US dollar [4, 8]. Most of the school/ 
institutional building in Nepal have shown lack of 
engineered design, lack of ductile detailing, 
structural deficiencies and low construction 
materials and techniques [4, 9-10]. These are the 
prime reasons which magnifies the importance to 
identify the performance level of institutional 
existing buildings that could be used in future 
explores to assure life safety during seismic 
events. 


In this context, this research focused on to study 
the behavior of institutional buildings under the 
effect of ground motions. The study identifies the 
status of damage grade: slight damage, moderate 
damage, extensive damage and complete damage 
of institutional building in synthetic earthquake. 
In fact, seismic risk analysis of existing 
institutional building is necessary to identify 
seismic vulnerability level under the effect of 
potential earthquake [11, 12]. It is achieved 
through numerical analysis. The building blocks 
are modelled in finite element-based software. 


The non-linear static and linear dynamic analyses 
are employed for numerical analysis. In dynamic 
analysis, building models are subjected to the 
synthetic accelerograms of the 2015 Gorkha 
earthquake. Based on the analyses, the analytical 
fragility curves are plotted in terms of probability 
of failure at every 0.1g interval of peak ground 
acceleration (PGA) with log normal distribution 
assumption. The results are highlighted for 


N. Ghimire and H. Chaulagain 


different seismic performance level in buildings 
for the earthquake of 475 years return period. 


2. Numerical Modelling of Structure: 


In this study, one of the school building of 
Pokhara University (School of Health and Allied 
Science, SHAS) was considered as a study 
institutional building. The building consists of 
three blocks. The building blocks model for 
numerical analysis was created and analyzed by 
ETABS software [13]. The software provided 
reliable result of linear and non-linear analysis 
that used to identify the behaviors of blocks under 
the static and dynamic loading. 


For numerical modelling, non-linear stress strain 
model for confined and unconfined concrete was 
considered as per Mander et al [14]. It consists of 
simple stress strain curve with controlling 
parameters for modelling of steel reinforcement. 
The beam and column sectional elements were 
modeled as elastic elements with plastic hinges at 
ends of members. The plastic hinges represent the 
concentrated behavior of structural member 
during numerical analysis. Default hinges 
characteristic used for the concrete sections was 
based on FEMA-356 [15] and ATC-40 [16] 
criteria. Flexural default hinge (M3) was assigned 
both ends of beams member and columns 
interacting (P-M2-M3) frame hinge type a 
coupled hinge property was assigned both lower 
and upper end of the member. 


2.1. Nonlinear Static Analysis 


In the recent years, non-linear static analysis 
gains the popularity among the researchers due to 
its simplicity and validity. Pushover analysis is 
conducted to estimate the strength and 
displacement capacity of existing structures [17]. 
Several international seismic guidelines proposed 
the standard procedures for pushover analysis 
[15, 16]. Generally, displacement controlled and 
force-controlled methods are used during 
pushover analysis. In the present study, 
displacement-controlled method was used for 
numerical analysis. During analysis, the building 
model was subjected to gravity load and 
monotonic displacement controlled lateral load 
patten in a seismic force continuously increasing 
from elastic to plastic behaviors of structure until 
the targeted displacement is exceeded. Target 
displacement of building is maximum 
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displacement of top floor level that was 
determined by standard guidelines [15-16]. The 
following steps were followed during analysis in: 


= create 3-D models of study institutional 
buildings 


= apply dead and live load on a numerical 
model 


= define plastic hinge properties on beam 
column element assuming 10% relative 
distance 


= apply the displacement on the top storey of 
structure whose value is large than those 
associated with target displacement 


= develop tables of roof displacement verse 
base shear or pushover curve 


2.2. Time History Analysis: 


The study of dynamic response of building 
structure behaviors under the specific loading 
(past earthquake) that may vary with time is 
known as time history analysis. The analysis was 
used to compute response of structure at a number 
of subsequent time instants. For the time history 
analysis, a representative Gorkha earthquake time 
history recorded at Kantipath used to evaluate 
structure response. Linear time history analysis 
compute the solution of dynamic equilibrium 
equation of building structure behaviors at given 
arbitrary time. The main purpose of linear time 
history analysis was to estimate displacement 
demand of building structure response in terms of 
top storey displacement. 
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Figure 1: Time history plot of 2015 Gorkha earthquake 
in Nepal 


To calculate structure demand parameter, it is 
necessary to determine the behavior of structures 
using the earthquake recorded data. In this study, 
synthetic earthquake data that consist of a 
simulated ground motion time history of Gorkha 
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earthquake was used. An appropriate set of 
acceleration was required for dynamic analysis of 
structure [18]. 


The earthquake data has compatible with elastic 
design response spectrum of this region [19]. 
First synthetic/artificial accelerograms were 
purposed by Gasparini and Vanmarke [20] which 
was based on design response spectrum. The time 
history data used in this study is presented given 
in Table 1, while the records are shown in Fig. 1. 


Table 1: Time history data used for dynamic 


analysis 
Name of Peak Ground 
Earthquake Acceleration (PGA) 
Gorkha earthquake 0.4g 


3. Case Study: 
3.1. Building Model: 


In this study, one of the school building of 
Pokhara University Nepal (School of Health and 
Allied Science, SHAS) was considered as a study 
building (see Fig. 2). The building is situated in 
Pokhara Metropolitan City, 30. The institutional 
building consists of three parts: left portion 


before the expansion joint is called block A, mid 
portion between two expansion joints is called 
block B and right portion after expansion joint is 
called block C (see Fig. 3-5). In this study, the 
response of the building can be observed by block 
wise. 


PoP rot 


Figure 2: (a) Existing SHAS building and (b) Plan of 
SHAS building 


Building block A, B and C measures 14.1 m 
x14.8 m, 24.6 m x 16.8 mand 14.1 m x 14.8 min 
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x and y direction, respectively. All the building 
blocks have storey height of 3.3 m and slab 
thickness 150 mm. Block A comprises five 
moment resisting frame in x and four moment 
resisting frame in y direction. Similarly, block B 
have six moment resisting frame in x and four 
moment resisting frame in y directions but in 
middle portion, one additional structural frame 


Table 2: Building section properties 


was inserted. Likewise, in building block C, five 
moment resisting frame in x direction and four 
moment resisting frame in y direction. The plan 
and 3D of the study building blocks structures are 
presented in Fig. 3-5. The details present of beam 
and column size and its reinforcement are in 
Table 2. 


Beam Column Slab 
Section Size (mm) Rebar Section Size (mm) Rebar 
3 O25 m 4 520 
B1 300x500 3 Ø20 C1- (G.F & 1* floor) 400x450 4016 
4 20 A 4 @20 
B2 300x500 3 Ø20 C2- (G.F & 1* floor) 400x450 8 Ø16 
5 20 af 4 O25 
B3 300x500 3 020 C3- (G.F & 1* floor) 400x450 4.920 
3 20 
B4 300x500 3920 C4- (G.F & 1* floor) 400x450 8 Ø 20 
B5 300x350 sie C5- (G.F & 1* floor) 400x450 8 O16 
4 @20 : 150 mm 
B6 300x350 i pap C1- (2™ & 3" floor) 400x450  8ø16 Thickness 
5 20 4 a 
B7 300x350 3920 C2- (2™ & 3 floor) 400x450 12 G16 
6 @20 ere 4 @25 
B8 300x350 3020 C3- (2™ & 3" floor) 400x450 4016 
5 @20 er 4 20 
B9 300x600 3 020 C4- (2™ & 3" floor) 400x450 4016 
6 20 nd rd 4 Ø16 
B10 300x600 3 Ø20 C5-(2"4 & 34 floor) 400x450 4012 


3.2. Material Properties: 


To identify behaviors of the structure, several 
material properties should be taken. The material 
properties of rebar and compressive strength of 
concrete were taken from experimental testing 
and blue print of existing drawing. The yield 
strength of rebar for all building blocks were 
adopted as 415 MPa. Schmidt hammer test is a 
non-destructive testing method that was 
performed to identify the compressive strength of 
concrete. Some specific material properties from 


experimental testing and blue print of existing 
drawing is tabulated in Table 3. During 
modelling, weight of infill wall was considered 
but stiffness of the masonry wall was ignored. 
Similarly, dead load was calculated as per 
drawing specified and live load was adopted as 
per IS-875 Part-II [21]. The live load in corridors, 
passages, staircase taken as 4 kN/m?. Similarly, 
live load in class room, store and bathroom taken 
as 3 kN/m?, 5 kN/m? and 2 kN/m’, respectively. 


Table 3: Concrete material properties of building. 


Properties 
Elastic modulus 
Shear modulus 
Weight density 
Poisson’s ratio 
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Values 
19364.92 MPa 
8068.72 MPa 

25 kN/m3 
0.2 
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(a) (b) 
Figure 3: Block A building model (a) Plan and (b) 3D view 
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(a) (b) 


(a) (b) 
Figure 5: Block C building model (a) Plan and (b) 3D view 
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4. Result and Discussion: 
4.1. Fragility Analysis: 


To derive the fragility curve from cumulative 
probability of failure correlate with increasing the 
demand displacement value (Sd) based on best 
fitted log-normal distribution function of 
equation that defined by median (Sc) and 
standard deviation parameters (£). The equation 
of probability of failure is defined as a function of 
Sd: 


(d) 


mene 


Pe) = 0 |=" 
Where, Æ ) stands for standard normal 
cumulative distribution function, Sd be the 
demand displacement and Sc and f be the 
medium of damage state value and standard 
deviation that represent total uncertainty. It is 


simply square root sum of the square’s 


combination of individual variability terms which 
is equivalent to 0.64 from HAZUS table [22]. 


e2 j 


(e) (f) 


Figure 6: Floor beam layout (a) Block A Ist floor, (b) Block B 1st floor, (c) Block C 1st floor, (d) Block A 2nd 
3rd 4th floors (e) Block B 2nd 3rd 4th floors (f) Block C 2nd 3rd 4th floors. 


(a) X- direction 


Base Shear (kN) 
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(a) 


Demand displacement was estimated by linear 
time history analysis of Gorkha earthquake. 
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(b) Y- direction 


Base Shear (KN) 


0 15 30 45 60 75 90 
Displacement (mm) 


(b) 
Figure 7: Pushover curve of Block A, B and C with (a) X-direction and (b) Y-direction 


Medium of damage states were estimated based 
on yield and ultimate displacement proposed by 
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Giovinazzi and Lagomarsino [12]. Idelization of 
capacity curve gives yield displacement (dy) and 
ultimate displacement (du) of building structure. 
The following damage states used to describe 
performance levels of building structure. 

=" Slight damage =0.7dy 
= Moderate damage = 1.5dy 
= Extensive damage = 0.5 ( dy + du) 


= Complete damage = du 


Table 4: Probability of failure of building 


Building Slight nen S Complete 
0 0 
(Block) (%) (%)  (%) (%) 
A 98 856 83 67.6 
B 97.9 802 774 60.2 
c 98.8 855 811 63.8 


Fragility curves are constructed with four damage 
states that represent probability of failure from 0 
to 1 corresponding to peak ground acceleration 
from 0 -1 g. Fig. 8 indicates the fragility curves 
of study building models which help 
engineer/decision maker to assure the seismic 
vulnerability condition of the structure and to 
predict the damage of building structure from 
possible future earthquake. Seismic hazard map 


A 
z 
= 
> 
= 
—] 
§ + Slight 
2 -=- Moderate 
3 Extensive 
a -Complete 
0 0.2 0.4 0.6 0.8 1 
PGA(g) 
(a) 


of Nepal proposed by Parajuli [23] indicates that 
PGA value of Pokhara for 10 % probability of 
exceedence in 50 years was expected to be 0.4 g. 
So, probability of failure of building at different 
damage state was observed at 0.4 g PGA. 


From Fig. 8, it can be seen that building block A 
has probability of failure in slight damage will be 
quite high, approximately 98 %. Whereas 
moderate, extensive and complete damage found 
as approximate 85.6 %, 83 % and 67.6 %, 
respectively. Similarly, building block B 
indicates that slight damage of building is higher 
as compared with other damage state. Slight, 
moderate, extensive and complete damage was 
estimated 97.9 %, 80.2 %, 77.4 % and 60.2 % 
respectively. In block C building model; slight, 
moderate, extensive and complete damage state 
was estimated 98.8 %, 85.5 %, 81.1 % and 63.8 
%, respectively. The above discussion result 
concluded that building performance levels 
depends upon the probability of failure of each 
damage state, if higher probability of failure of 
damage states less performance levels of 
buildings. Finally, more than 50 % damage state 
was estimated of all building block are present in 
Table 4. 


-+ Slight 
-® Moderate 
~*- Extensive 


—* Slight 


-= Moderate 
-+ Extensive 
~~ Complete 


-Complete 


(b) (c) 


Figure 8: Fragility curves (a) block A, (b) block B and (c) block C 


5. Conclusions: 


The study focused to construct the fragility curve 
of RC institution building. For this, RC 
institutional building in Pokhara University 
(SHAS) was considered as case study. For 
obtaining fragility curve linear and non-linear 
analyses were performed with finite element 
software. Fragility curves were plotted with the 
probability of failure at every 0.1 g interval of 
PGA that used to estimate building structure 
damage/performance levels under the effect of 
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potential earthquake. The conclusions of 
analytical fragility curve can be summarized as: 


= Seismic fragility curve of all the building 
models are constructed with four damage 
state (slight, moderate, extensive and 
complete). An increase in PGA leads to an 
increased fragility. 


=" Fragility curve of building depends upon 
building properties and configuration but 
not an earthquake data. 
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= The slope of the fragility curve was large at 
lower PGA and smaller at higher PGA. 
Finally, more than 50 % of all damage states 
were found at 0.4 g PGA value which 


concluded 


that a building was low 


performance level at 475 years return period 
of earthquake. 
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